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Abstract 

Given a family of k disjoint connected polygonal sites in general position and of total complexity n, 
we consider the farthest-site Voronoi diagram of these sites, where the distance to a site is the distance 
to a closest point on it. We show that the complexity of this diagram is 0(n), and give an 0(n log^ n) 
time algorithm to compute it. We also prove a number of structural properties of this diagram. In 
particular, a Voronoi region may consist of fc — 1 connected components, but if one component is 
bounded, then it is equal to the entire region. 

1 Introduction 

Consider a family S of geometric objects, called sites, in the plane. The farthest-site Voronoi diagram 
of S subdivides the plane into regions, each region associated "with one site P € S, and containing those 
points x G for "which P is the farthest among the sites of S. 

While closest-site Voronoi diagrams have been studied extensively |1] , their farthest-site cousins have 
received some"what less attention. For the case of (possibly intersecting) line segment sites, Aurenham- 
mer et al. [3] recently presented an O(nlogn) time algorithm to compute their farthest-site diagram. 

Farthest-site Voronoi diagrams have a number of important applications. Perhaps the most "well- 
kno"wn one is the problem of finding a smallest disk that intersects all the sites [T]. This disk can be 
computed in linear time once the diagram is kno"wn, since its center is a vertex or lies on an edge of the 
diagram. Another standard application is to build a data structure to quickly report the site farthest 
from a given query point. 

We are here interested in the case of complex sites "with non-constant description complexity. This 
setting "was perhaps first considered by Abellanas et al. [Tj: their sites are finite point sets, and so the 
distance to a site is the distance to the nearest point of that site. Put differently, they consider n points 
colored "with k different colors, and their farthest- color Voronoi diagram subdivides the plane depending 
on "which color is farthest a"way. The motivation for this problem is the one mentioned above, namely to 
find a smallest disk that contains a point of each color — this is a facility location problem "where the goal 
is to find a position that is as close as possible to each of k different types of facilities (such as schools, post 
offices, supermarkets, etc.). In a companion paper [5] the authors study other color-spanning objects. 

The farthest-color Voronoi diagram is easily seen to be the projection of the upper envelope of the 
k Voronoi surfaces corresponding to the k color classes. Huttenlocher et al. [TT] sho"w that this upper 
envelope has complexity &{nk) for n points, and can be computed in time 0{nk logn) (see also the book 
by Sharir and Agarwal [THl §8.7]. 
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Figure 1: (a) The bisector of two polygons can be a closed curve, (b) The medial axis 9H(Pi) outside 
of Pi. (c) The Voronoi region TZ{Pi), shown in light gray, (d) The farthest-polygon Voronoi diagram 

mPi,P2})- 

Van Kreveld and Schlechter [T^ consider the farthest-site Voronoi diagram for a family of disjoint 
simple polygons. Again, they are interested in finding the center of the smallest disk intersecting or 
touching all polygons, which they then apply to the cartographic problem of labeling groups of islands. 
Their algorithm is based on the claim that this farthest-polygon Voronoi diagram is an instance of the 
abstract farthest-site Voronoi diagram defined by Mehlhorn et al. [15| — but this claim is false, since the 
bisector of two disjoint simple polygons can be a closed curve, see Fig. [TJa). In particular, Voronoi 
regions can be bounded (which is impossible for regions in abstract farthest-site Voronoi diagrams). 

Note that the farthest-polygon Voronoi diagram can again be expressed as the upper envelope of k 
Voronoi surfaces — but this does not seem to lead to anything stronger than near-quadratic complexity 
and time bounds. 

We show in this paper that, in fact, the complexity of the farthest-polygon Voronoi diagram of k 
disjoint simple polygons of total complexity n is 0{n). We also show some structural properties of this 
diagram. In particular, Voronoi regions can be disconnected, and in fact, the region of a polygon P can 
consist of up to fc — 1 connected components. However, if one connected component is bounded then 
it is equal to the entire region of P; moreover, the region is simply connected and the convex hull of P 
contains another polygon in its interior. Furthermore, the Voronoi regions consist, in total, of at most 
2k — 2 connected components, and this bound is tight. 

Algorithms for computing closest-site Voronoi diagrams make use of the fact that Voronoi regions 
surround and are close to their sites. Similarly, algorithms for computing farthest-site Voronoi diagrams 
make use of the unboundedness of the Voronoi regions, and often build up regions from infinity [3]. 
The difficulty in computing farthest-polygon Voronoi diagrams is that neither of these properties holds: 
Voronoi regions can be boimded, and finding the location of these bounded regions is the bottleneck in 
the computation. 

We give a divide-and-conquer algorithm with running time 0{n\og^ n) to compute the farthest- 
polygon Voronoi diagram. Our key idea is to build point location data structures for the partial diagrams 
already computed, and to use parametric search on these data structures to find suitable starting vertices 
for the merging step. This idea may find applications in the computation of other complicated Voronoi 
diagrams. Our algorithm implies an 0{nlog^ n) algorithm to compute the smallest disk touching or 
intersecting all the input polygons. 

We note that for a family of disjoint convex polygons, finding the smallest disk touching all of them 
is much easier, and can be solved in time 0{n), where n is the total complexity of the polygons [12] . 

In Section [21 we start with some preliminaries and give a definition of farthest-polygon Voronoi 
diagrams. We prove, in Section [31 some properties on the structure of these diagrams and bound their 
complexity. In Scction[31 we present an algorithm for computing such diagrams and conclude in Section[3 
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Figure 2: A polygonal site P with three bounded pockets and seven unbounded ones, and its convex 
hull. A point x and the disk Dp{x). 

2 Preliminaries 

We consider a family S oi k pairwise-disjoint polygonal sites of total complexity n. Here, a polygonal 
site of complexity m is the union of m line segments, whose relative interiors are pairwise disjoint, but 
whose union is connected, see Fig. [5) (In other words, the corner^ and edges of a polygonal site form 
a one-dimensional connected simplicial complex in the plane.) In particular, the boundary of a simple 
polygon is a polygonal site. For a point x G R^, the distance d{x, P) between x and a site P S 5 is the 
Euclidean distance from x to the closest point on P. 

A pocket P of P is a connected component of CII(P) \ P, where CII(P) is the convex hull of P. A 
pocket V is hounded if it coincides with a bounded connected component of \P, unbounded otherwise. 

The features of a site P are its corners and edges. We say that a disk touches a corner if the corner 
lies on its boundary. A disk touches an edge if the closed edge touches the disk in one point and if its 
supporting line is tangent to the disk. A disk touches a site P if the disk touches some of P's features 
and if the disk's interior does not intersect P. 

We assume that the family S is in general position, that is, no disk touches four features, no line 
contains three corners, and no two edges are parallel. 

Before defining the farthest-polygon Voronoi diagram of a family of polygonal sites, we define the 
medial axis of a polygonal site. 

Medial axes. For a site P G 5, we define the function *p : ^ r as *p(.t) = d{x,P). The 
graph of ^'p is a Voronoi surface] it is the lower envelope of circular cones for each corner of P and of 
rectangular wedges for each edge of P. The orthogonal projection of this surface on the plane induces a 
subdivision of the plane: each 2D cell of this subdivision corresponds to a feature w of P, and it is the 
set of all points a; G such that w is or contains the unique closest point on P to x (here, edges of P 
are considered relatively open, so the cell of a corner is disjoint from the cells of its incident edges). The 
medial axis of P, denoted S[Jl(P), consists of the arcs and vertices formed by the boundaries between 
these cells. By extension, we call the cells of the subdivision the cells of the medial axis. 

For a point x G and a site P, let Dp{x) denote the largest disk centered at x whose interior does 
not intersect P (and which is therefore touching P, see Fig. [5]). If Dp{x) touches P in a single feature 
w, then X lies in the cell of the medial axis subdivision associated with w. If Dp{x) touches P in two 
different features, then x lies on an arc of 9Jl(P), and if Dp{x) touches P in three different features, then 
x is a vertex of 9Jl(P). Note that 9Jl(P) contains some special arcs, called spokes, that separate the cell 
of a corner from the cell of an incident edge (sec Fig.[TJb)); a spoke arc is the locus of centers of circles 
that intersect P only at a corner and that are tangent to the line supporting one of its incident edgesH 

The medial axis 93t(P) restricted to \ P forms a forest. By this definition, the arc cndpoints that 
lie on P (at a corner) are not part of the forest; we consider nonetheless these cndpoints to be leaves 
of the forest. It follows that several leaves may coincide at a corner. More precisely, each convex angle 

^We reserve the word vertex for vertices of the Voronoi diagram. 

■^Note that the medial axis of P is often defined as the locus of the centers of circles that are tangent to P in two or 
more points. Such a definition would not include spokes as part of the medial axis. 
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around a corner induces a leaf at that corner, and each reflex angle around a corner induces two leaves 
incident to two spokes at that corner. Spokes are always incident to a leaf at a corner. Notice that the 
leaves always lie at the corners of P or at infinity. 

The medial axis 9Jt(P) consists of exactly one tree for each pocket of P, and two isolated spokes (with 
endpoints at infinity) for each edge of P appearing on CH(P). The tree T of an unbounded pocket P 
contains exactly one infinite arc, all other leaves of T are corners of V. The leaves of the tree T of a 
bounded pocket V are exactly the corners of V. 

Farthest-polygon Voronoi diagrams. We now consider the function $ : n> E defined as = 
nia.xp^s p{x). The graph of $ is the upper envelope of the surfaces for P E S. The surface 
<f> consists of conical and planar patches from the Voronoi surfaces and the arcs separating such 
patches are either arcs of a Voronoi surface p (we call these medial axis arcs) , or intersection curves of 
two Voronoi surfaces 'I'p and (we call these pure arcs). These arcs are hyperbolic arcs that lie in 
vertical planes, parabolic arcs, or straight-line segments. They correspond respectively to the intersection 
of two cones, a cone and a plane, and two planes. The vertices of $ are of one of the following three 
types: 

• Vertices of one Voronoi surface ^! p. We call these medial axis vertices. 

• Intersections of an arc of ^'p with a patch of another surface We call these mixed vertices. 

• Intersections of patches of three Voronoi surfaces ^'p, 'I'p. We call these pure vertices. 

The projection of the graph of onto the plane induces the farthest-polygon Voronoi diagram ^{S) of 
S. It is a subdivision of the plane into cells, arcs, and vertices. The arcs are either parabolic or straight, 
since hyperbolic arcs that lie in vertical planes project into line segments. Each cell corresponds to a 
feature w of a site P, the feature is the nearest among the features of P, but is further away than the 
nearest feature of any other site. The arcs and vertices of d{S) are the orthogonal projections of the 
arcs and vertices of $ and they inherit their types (pure, mixed, and medial axis). The farthest-polygon 
Voronoi diagram is therefore completely analogous to the farthest-color Voronoi diagram [T] . 

For a point x £ R^, let D{x) = Ds{x) denote the smallest disk centered at x that intersects all sites 
P e 5. By definition, there is always at least one site that touches D{x) without intersecting its interior, 
and the radius of D(x) is equal to <I>(x). By our general position assumption, only the following five 
cases can occur. 

• If D{x) touches one site P in only one feature w, and all other sites intersect the interior of D{x), 
then X lies in a cell of d{S), and the cell belongs to the feature w of P. 

• If D{x) touches one site P in two or three features, and all other sites intersect the interior of D{x), 
then X lies on a medial axis arc or medial axis vertex of 5('5), and is incident to cells belonging to 
different features of P. 

• If D{x) touches one feature w of site P, one feature u of site Q, and all other sites intersect the 
interior of D{x), then x lies on a pure arc separating cells belonging to features w and u. 

• If -D(x) touches two features of site P and one feature of site Q, and all other sites intersect the 
interior of D{x), then x is a mixed vertex incident to a medial axis arc of P. 

• If D{x) touches one feature each of three sites P, Q, and R, and all other sites intersect the interior 
of D{x), then a; is a pure vertex. 

Put differently, vertices of ^{S) are points x where D{x) touches three distinct features of sites. If 
all three features are on the same site, the vertex is a medial axis vertex. If the three features are on 
three distinct sites, then the vertex is a pure vertex. In the remaining case, if two features are on a site 
P, and the third feature is on a different site Q, the vertex is a mixed vertex. 

Fig. [3] illustrates all different vertex types. Figs.[3Ua) and[3Jb) show the possible medial axis vertices; 
they differ in whether the triangle formed by the three features contains the vertex or not. Similarly, 
Figs.jSl^c) and[31^d) show the pure vertex types. The bottom row shows the three possible types of mixed 
vertices. Again, we have to distinguish whether the three features enclose the vertex or not, and in the 
latter case we need to distinguish which two features are on the same site. 

Consider an arc a of ^{S). If a point x moves continuously along a, then <^{x) — which is the radius 
of D{x) — changes continuously. The local shape of ^{S) in a neighborhood of a vertex v is determined 
solely by the features defining the vertex. For each type of vertex shown in Fig. [21 we can therefore 
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Figure 3: The different types of vertices in the farthest-polygon Voronoi diagram. Pure arcs are shown 
solid, medial axis arcs dashed. The arrows indicate the direction of increasing $(a;) in a neighborhood 
of the vertex. 

uniquely determine whether 4>(a;) increases or decreases along an arc in a neighborhood of the vertex. 
We indicate the direction of increasing <I>(x) along an arc by arrows in the figure. 

In addition to the seven types of vertices discussed above, we need to consider vertices at infinity, 
that is, we consider the semi-infinite arcs of ^{S) to have a degree-one vertex at their end. For a vertex 
V at infinity, the "disk" D(v) is a halfplane, and we have two cases: 

• If D{v) touches one site P in two features, and all other sites intersect the interior of D{v), then 
D{v) is the "infinite" endpoint of a medial axis arc, and we consider it a medial axis vertex at 
infinity. 

• If D{v) touches two distinct sites, and all other sites intersect its interior, then D{v) is the "infinite" 
endpoint of a pure edge, and we consider it a pure vertex at infinity. 

Finally, we define the Voronoi region of a site P £ S. The Voronoi region TZ{P) of P is simply 
the union of all cells, medial axis arcs, and medial axis vertices of belonging to features of P. 

Voronoi regions are not necessarily connected, as we will see in the next section. We call each connected 
component of a Voronoi region a Voronoi component. 

3 Structure and complexity 

In this section, we prove some properties on the structure of farthest-polygon Voronoi diagrams of 
polygonal sites and we bound their complexity. Farthest-polygon Voronoi diagrams contain three different 
types of vertices, as defined in the previous section. Note first that the number of medial axis vertices is 
bounded by the total complexity of all medial axes, which is 0{n). 

Now, we first bound the number of mixed vertices. For that purpose, we show that when a tree T 
of 9Jl(P) intersects the Voronoi region TZ{P) oi P G S, then the intersection T n TZ{P) is a connected 
subtree. We start with a preliminary lemma. 

Lemma 1. Let ^ he a path in 9Jt(P), let Q £ S \ {P} be another site, and let 7q be the part of "/ that 
is closer to Q than to P. Then 7q is a connected subset of 7, that is, a subpath. 

Proof. We can assume 7 to be a maximal path in T, connecting a corner w of P with another corner u 
or a medial-axis vertex at infinity. Assume for a contradiction that there are points x, y, 2; on 7 in this 
order such that x, z g 7q, but y ^ ^q. 
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We first consider the case where y lies on a spoke of 9Jt(P) induced by a corner c of P (and one of 
its incident edges). The spoke is incident to a leaf of 9JT(P) located at c and, without loss of generality, 
X lies on this spoke between c and y. Since y ^ Jq, the disk Dp{y) centered at y and touching c does 
not intersect Q. The disk Dp{x) is included in Dp{y) and thus, it does not intersect Q either, which 
contradicts our hypothesis. 

We now consider the case where y docs not lie on a spoke of Let 51 be the connected component 

of \ P that contains y. Since y does not lie on a spoke of 9Jl(P), the disk Dp{y) touches P in fc 2 
distinct points (fc ^ 3 by the general position assumption), and thus Dp{y) partitions U, into fc + 1 
connected components: Dp{y) and k other components denoted Ai, . . . , A^- Since P is a polygon, the 
structure of its medial axis is well understood. In particular, in the neighborhood of y, 3Jl(P) consists 
of k arcs (straight or parabolic), and for every point p on any single one of these arcs, Dp{p) intersects 
one and the same components Ai, and is contained in AiU Dp{y). Also, since ?OT(P) consists of k arcs in 
the neighborhood of y, point y splits the medial axis tree T that contains 7 into k subtrees Ti, ■ ■ ■ ,Tk- 

Now, we observe that any open disk D, that does not intersect P, cannot contain points in two 
distinct components Ai and Aj. Indeed, the boundary of D would have to intersect the boundary of 
Dp(y) in at least four points, implying that the two disk coincide, and thus that D intersects neither Ai 
nor Aj (since D is open). 

For any p £ 7i, distinct from y, Dp{p) is included in Q but not in Dp{y). Thus, the interior of Dp{p) 
intersects [Jj Aj but not P. Hence, it intersects only one of the Aj. Furthermore, by continuity, the 
interior of Dp{p) intersects the same Aj for all p in Ti- We can thus assume without loss of generality 
that, for all p in 71, the interior of Dp{p) intersects Ai and none of the other Aj. Since Dp{p) lies in fi, 
it also follows that, for all p in 71, Dp{p) hes in Ai U Dp{y). 

Now, X and z belong to two distinct subtrees, say 71 and 72, respectively. Thus, Dp(x) lies in 
Ax U Dp{y) and Dp{z) lies in A2 U Dp{y). By assumption, both Dp{x) and Dp{z) intersect Q, thus Q 
intersects both Ai Li Dp{y) and A2 U Dp{y). Hence, since Q is connected and does not intersect P, Q 
must intersect Dp(y), which is a contradiction and concludes the proof. □ 

Lemma 2. Let T be a tree o/9}l(P). Then T r\TZ{P) is a connected subtree of T ■ 

Proof. Pick two points p,q € T H TZ{P), let 7 be the path on T from p to q, and let x be any point 
between p and q on 7. We need to show that x G TZ{P), which is equivalent to showing that x is closer 
to every site Q P than to P. Let Q be such a site. Since p,q G TZ{P), p and q are closer to Q than 
to P. By Lemma [2 this implies that x is closer to Q than to P. □ 

The lemma allows us to bound the number of mixed vertices of S{S). 

Lemma 3. The number of mixed vertices of ^(S) for a family of disjoint polygonal sites of total com- 
plexity n is 0{n). 

Proof. Consider a site P S 5 of complexity m. Its medial axis 9Jl(P) has complexity 0{m). By LemmajU 
for each tree T of S!7l(P) the intersection T n TZ{P) is a connected subtree. Since the mixed vertices 
on T are exactly the finite leaves of this subtree, this implies that the number of mixed vertices on 9Jl(P) 
is 0(m). Summing over all P G iS then proves that the number of mixed vertices of 5^(5) is 0(n). □ 

We next consider the vertices at infinity. 

Lemma 4. The number of pure vertices at infinity of ^{S) is at most 2k — 2. The total number of 
vertices at infinity of^{S) is 0(n). 

Proof. For two sites P,Q<eS, consider the diagram d{{P, Q})- A pure vertex at infinity corresponds to 
an edge of CH(P U Q) supported by a corner of P and a corner of Q. But CH(P U Q) can have at most 
two such edges, since P and Q are disjoint and both are connected, and so d{{P, Q}) has at most two 
pure vertices at infinity. 

Consider now again 3^(5), and let cr(S) denote the sequence of sites whose Voronoi regions appear at 
infinity in circular order, starting and ending at the same region. We claim that a{S) is a Davenport- 
Schinzcl sequence of order 2, and has therefore length at most 2fc — 1 [T^. Indeed, cr(iS) has by definition 
no two consecutive identical symbols. Assume now that there are two sites P and Q such that the 
subsequence PQPQ appears in (j{S). If we delete all other sites, then <t({P, Q}) would still need to 
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contain the subsequence PQPQ, and therefore d{{P, Q}) would contam at least three pure vertices at 
infinity, a contradiction to the observation above. 

It now suffices to observe that the pure vertices at infinity are exactly the transitions between con- 
secutive Voronoi regions, and their number is at most 2k ~ 2. 

All remaining vertices at infinity are medial axis vertices. Since the total complexity of all medial 
axes is 0{n), the bound follows □ 

We proved so far that the number of mixed and medial axis vertices is 0{n) and, furthermore, that 
there are at most 2k — 2 pure vertices at infinity. It remains to bound the other pure vertices, for which 
we first need to prove a few basic properties. 

We start by discussing a monotonicity property of cells of ^{S) . Let C be a cell of 5^(5) belonging to 
feature w of site P. For a point x <E C, let x* be t he ^ oint on w closest to x. Let fx be the directed line 
segment starting at x and extending in direction x*x until it reaches 97l(P) (a semi-infinite segment if 
this does not happen). We call fx the fiber of x. We note that if w is an edge, then all fibers of C are 
parallel, and normal to w; if it; is a corner then all fibers are supported by lines through w. 

Lemma 5. For any x G C, the fiber fx lies entirely in C (and therefore in TZ{P) )■ 

Proof. The disk D{x) touches P in x* only, and its interior intersects all other sites. When wc move a 
point y from x along fx, the disk D centered at y through x* keeps containing D{x), and it therefore 
still intersects all other sites. This implies that y S C as long as D does not intersect P in another point. 
This does not happen until we reach 3Jl(P). □ 

An immediate consequence, which we will use for computing Voronoi diagrams (Section U]), is that 
cells are "monotone": 

Lemma 6. The boundary of a cell C of ^{S) belonging to feature w consists of two chains monotone 
with respect to w, that is, monotone in the direction of w if w is an edge, and rotationally monotone 
around w if w is a corner. The lower chain is closer to the feature and consists of pure arcs only, the 
upper chain consists of medial axis arcs only. 

Proof. Let P be the site containing the feature w. Consider a half-line £ with origin on w, and normal 
to ui if u) is an edge. Let x be the point closest to w in £ DC. It is straightforward that the entire fiber 
fx lies in C, and no point z on £ beyond the medial axis can be in C since the feature of P closest to 
any such z cannot be w. Hence, the boundary of C consists of two monotone chains with respect to w. 
Moreover, the upper chain consists of medial axis arcs by definition of the fibers fx- The lower chain 
consists of pure arcs, because if a point x on the lower chain was on the medial axis, then the fiber fx 
would be reduced to point x, by definition; thus x would also be on the upper chain, implying that x is 
an endpoint of the two chains. □ 

We now show (Lemma [5]) that if a Voronoi region is bounded, then it is connected (we actually show 
that it is simply connected, but will not use that fact in this paper). This property is tight in the sense 
that, as shown in Lemma 1 101 a single Voronoi region may consist of up to fc — 1 unbounded connected 
components; we postpone the proof of this property to the end of the section. 

Lemma 7. // a connected component of the Voronoi region 'R-{P) of a site P E S is bounded, then P 
properly contains another site inside one of its pockets. 

Proof. Let C be a bounded connected component oi'R{P). We first observe that C contains some points 
of the medial axis 2Tt(P) of P. Indeed, let a; G C and consider its fiber fx- By Lemma [S] and since C is 
bounded, fx does not extend to infinity and, therefore, one of its endpoints lies on 9Jl(P). 

Let a; be a point in C n S[)T(P). If x lies in a pocket of P that does not share an edge with CII(P) 
(that is, the pocket is a hole in P), then this pocket docs contain all the other sites in 5 \ {P} and the 
lemma is proven. Otherwise, we let x move along the medial axis 93t(P) up to infinity. At some point x' , 
the point must exit from the bounded region C. This means that Ds{x') = Dp(x') is tangent to another 
site Q and does no longer intersect it properly. It then follows from the fact that x' lies on the medial 
axis of P that the site Q lies entirely in the pocket of P associated with the tree of 07l(P) containing x 
and x' . □ 
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Lemma 8. // a connected component of the Voronoi region TZ{P) of a site P € 5 is bounded, then TZ{P) 
is simply connected. 

Proof. By Lemma [71 P contains another site Q in one of its pockets *p. Let T be the tree of 2H(P) that 
corresponds to *p. 

Let X be any point in TZ{P). The disk D(x) touches P in a point x*, and its interior intersects all 
other sites, including Q. Hence, D{x^ properly intersects the pocket Cp. It follows that when moving a 
point y from x in the direction of x*x, the disk centered at y through x* keeps containing D{x), and it 
therefore keeps intersecting *p. This implies that, at some finite point y = x', the disk D centered at 
x' and tangent to P at x* becomes tangent to P at some other point, hence x' lies on 9n(P). Moreover, 
x' lies on the tree T because the disk D properly intersects *p. 

Hence, for any point x in TZ{P), the fiber fx is a segment joining x to a point x' on T. Fmthermore, 
the fiber fx lies in TZ(P), by Lemma [5l and Tr\TZ{P) is a connected tree, by Lemma [5] Therefore, TZ{P) 
is connected. 

It remains to show that Tl{P) is simply connected. We have shown that for any x G TZ{P), the fiber 
fx is a segment contained in TZ{P), and connecting x to x' € T H 7?,(P). By moving the points of TZ{P) 
along their fiber, we can design, as follows, a continuous deformation retraction of 7^(P) onto the tree 
T n TZ{P) which implies that Ti-iP) is simply connected. 

More precisely, it easy to check that the map F : TZ{P) x [0,1] — >■ R^, {x,t) h- > {1 — t)x + tx' is 
continuous. Furthermore, P is a (strong) deformation retraction since, for all x e TZ{P), t e [0,1] and 
meTn 7^(P), we have F{x, 0) = x, F{x, 1) e T n 7^(P), and P(m, t) = m. We have thus exhibited a 
deformation retraction of TZ{P) onto Tn TZ{P), which implies that these two point sets have the same 
fundamental group [lOi Proposition 1.17]. Since TnTZ{P) is a tree, its fundamental group is trivial and 
so is the fundamental group of TZ{P). Any pair of points x and y in TZ{P) can be connected with a path 
made of their fibers fx and fy and a path in Tn TZ{P), thus 7^(P) is path connected. Together with 
having a trivial fundamental group, this property makes TZ{P) simply connected |10l p. 28]. □ 

We can now conclude our analysis of the complexity of farthest-polygon Voronoi diagrams. 

Theorem 9. The farthest-polygon Voronoi diagram of a family of k disjoint polygonal sites of total 
complexity n has 0{k) pure vertices and total complexity 0{n). It consists of at most 2k — 2 Voronoi 
components and this bound is tight in the worst case. 

Proof. The farthest-polygon Voronoi diagram contains three different kinds of vertices. The number of 
medial axis vertices is clearly only 0(n), since the total complexity of all 9H(P) for P G 5 is only 0{n). 
In Lemma [21 we showed that the number of mixed vertices is also only 0{n). It remains to bound the 
number of pure vertices of ^{S). 

Let ki be the number of bounded Voronoi components. By Lemma [SI each of these components 
corresponds to a different site and only the remaining k — fci sites can contribute to form vertices at 
infinity. By the proof of Lemma [H there are at most 2(fc — fci) — 2 pure vertices at infinity, and therefore 
at most 2{k — ki) —2 unbounded Voronoi components. It follows that the total number of Voronoi 
components is at most 2/c — fci — 2 ^ 2fc — 2. Moreover, the construction of Lemma [TUl shows that this 
bound is tight. 

Let us now consider the graph G formed by the pure arcs and pure vertices of 5^(5) . Mixed vertices 
appear as vertices of degree two in G (see Fig. [3]), medial axis vertices do not appear at all. The faces 
of G are exactly the Voronoi components. Since G has at most 2fc — 2 faces, Euler's formula implies that 
G has 0{k) vertices of degree three. □ 

Finally, we prove that, as mentioned above, a single Voronoi region of S{S) can have up to A; — 1 
connected components. 

Lemma 10. A single Voronoi region of ^(S) can have fc — 1 connected components and this bound is 
tight. 

Proof. The construction is shown in Fig.[31 It consists of one (fc — l)-regular polygon R and k — 1 polygonal 
chains Ci , . . . , Ck-i. Let ei , 62, . . . , Ck-i denote the edges of R in circular order. We inductively construct 
the polygonal sites Gi, i = 1,2, . . . ,k — 1 as follows. For the supporting line I of Ci, let be the closed 
halfplane containing R and l~ be the other. Then, consider the intersection G* between l~^ and the four 
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Figure 4: The Voronoi region of polygon R has k — I connected components. 



edges of a square that contains i?, Ci, . . . , Ci-i inside. We define C; as the set of points of C* whose 
distance to I is larger than some fixed small e > 0. Ci has at most four edges. Note that contains 
Ci completely. Consider a ray from to infinity in l~ which is orthogonal to I. Since l~ intersects all 
sites but Ci, the endpoint at infinity of this ray lies in the region Ci — TZ{Ci). On the other hand, for a 
sufficiently small e, there is a line L passing through Vi (the vertex incident to e^^i and e^) such that we 
can define as an open halfplane containing i?\{7Ji} and L~ as the other open halfplane intersecting 
all the other sites but R. The endpoint at infinity of the ray from Vi to infinity in L~ which is orthogonal 
to L lies in a connected component of TZ{R), which we call TZi. 

Therefore at infinity 7?.i, Ci, 7^2, C2, • • • ,7^fc_i,Cfe_i appear in turn. Note that for a point x in the 
region TZ{R), its fiber is an infinite ray because R is convex. For i = 1,2, . . . ,k — 1, consider the half- 
line Xi(j)i from a point at infinity (jji G C; to a point Xi G R closest to cjji. If x G Xi(f)i lies in TZ{R), then 
fx C TZ{R), which is impossible because (j>i G Ci. Define W = RLl[J^ Xi(f)i. We then have W n Ti[R) = 0, 
and \ consists of — 1 unbounded connected subsets of the plane. The connected subset bounded 
by 01+1 and R contains TZi completely. It follows that TZi ^ TZj when i j. 

It remains to show that the bound of A: — 1 is tight. By Lemma SI there are at most 2k — 2 pure 
vertices at infinity, and thus at most 2k — 2 unbounded Voronoi components. Hence, one single Voronoi 
region has at most k — I unbounded Voronoi components (since two neighboring components cannot 
correspond to the same site). This concludes the proof because, by Lemma |51 if a Voronoi region has a 
bounded component, then the entire region is connected. □ 

4 Algorithm 

The proof of Theorem [5] suggests an algorithm for computing the Voronoi diagram by sweeping the arcs 
of the graph G. This is roughly equivalent to computing the surface $ by sweeping a horizontal plane 
downwards, and maintaining the part of $ above this plane. This is essentially the approach used by 
Aurenhammer et al. [3] for the computation of farthest-segment Voronoi diagrams. However, this does 
not seem to work for our diagram because of the mixed vertices of type (f) (see Fig. [3]), where $ has a 
local maximum. We instead offer a divide-and-conquer algorithm. 

Theorem 11. The farthest-polygon Voronoi diagram ^{S) of a family S of disjoint polygonal sites of 
total complexity n can be computed in time 0(nlog"^n). 

Proof. Let S = {Pi, . . . ,Pk}, and let rii be the complexity of P^. If fc = 1, then ^{S) is simply the medial 
axis Vyi{Pi), which can be computed in time O(nlogn) [2]. Otherwise, we split S into two disjoint 
families iSi,52 as follows: 

• If there is a site Pi with complexity Ui ^ n/2, then Si = {Pi} and S2 = S \ {Pi}. 

• Otherwise there must be an index j such that ?i/4 ^ J2i=i ^ 3?2/4. We let Si = {Pi, . . . , Pj} 
and 52 = {P,+i, • ■ • ,Pfe}. 
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d{Si) d{S2) Purple chains ^{S) 

Figure 5: Merging ^{Si) and 3^(52) to obtain 5(5). 



Wc recursively compute S'(iSi) and ^{S2)- We show in the rest of this section (see Lemma [TT]) that we 
can then merge these two diagrams to obtain ^(S) in time 0(nlog^ ri), proving the theorem. □ 

We now discuss the merging step. We are given the farthest-polygon Voronoi diagrams i?(5i) and 
^{82) of two families Si and ^2 of pairwise disjoint polygonal sites, and let S = iSi U 52- 

Consider the diagram S{S) to be computed. We color the Voronoi regions of 5(5) defined by sites in 
Si red, and the Voronoi regions defined by sites in ^2 blue. A pure arc of 5(5) is red if it separates two 
red regions, and blue if it separates two blue regions. The remaining pure arcs, which separate a red and 
a blue region, are called purple. A vertex of 5(5) is purple if it is incident to a purple arc. We observe 
(see Fig. [3]) that, by our general position assumption, every purple vertex not at infinity is incident to 
exactly two purple arcs, and so the purple arcs form a collection of bounded and unbounded chains (see 
Fig.©. 

As we will see in Lemma 1171 merging 5(5i) and 5(52) can be done in linear time once all purple 
arcs are known, because the diagram 5(5) consists of those portions of 5(5i) lying in the red regions of 
5(5), and those portions of 5(52) lying in the blue regions of 5(5), see Fig. [SJ 

We show below how the purple chains of 5(5) can be computed in time 0(n log^ n). We first show 
how to compute at least one point on every chain and then how to "trace" a chain from a starting point. 

Lemma 12. The vertices at infinity of the purple chains can be computed in time O(nlogn). 

Proof. We show how to compute all the pure vertices at infinity in time O(nlogri). There are 0{k) such 
pure vertices (by Lemma 2]), and we can easily deduce the purple ones from those in 0{k) time. 

For site G 5 and angle (j) G [0, 2tt), let £i{(j>) be the oriented line with direction (f) tangent to Pi and 



keeping Pi entirely on its left, see Fig. 6(a) Let gi{(l)) be the signed distance from the origin to £i{(j)) 
(positive if the origin lies left of £i{4>), negative otherwise). If Pi is a polygonal site of complexity m, we 
first compute the convex hull CH(Pi) in time O(TOlogTO), and we deduce a description of the function 
gi in time 0{m). 

We then compute the lower envelope g of the functions gi. The pure vertices at infinity correspond 
exactly to the breakpoints of this lower envelope, since they correspond to half-planes (or disks with 
centers at infinity) touching two sites and whose interiors intersect all other sites. Such a half-plane is 



illustrated in gray in Fig. 6(a) Moreover, two functions gi, gj can intersect at most twice since each 
intersection corresponds to a pure vertex at infinity of 5({-Pi, Pj}) which admits at most two such vertex as 
argued in the proof of Lemma|31 Hence, the lower envelope can be computed in time 0(nlog?T,) [19] . □ 



Lemma 13. Any bounded purple chain contains a mixed vertex o/5(5). 

Proof. A bounded purple chain is a compact set in the plane, and so the restriction of $(•) to such a 
chain admits a maximum. Such a maximum appears at a vertex, denoted v. Indeed, since an arc a of 
5(5) is defined by two features (corners or edges), the graph of $(•) restricted to a is the intersection 
of the two Voronoi surfaces — which are cones or wedges — induced by the two features; thus $(•) cannot 
have a local maximum in the interior of a (it may have a local minimum). 

Consider now the arcs of 5(5) oriented, in a neighborhood of their endpoints, in the direction of 
increasing $(•). Then, the two purple arcs incident to v point toward v. Now observe that purple arcs 
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4(0) =4(0) 4(0) 4(0) 




Figure 6: (a) Constructing the pure vertices at infinity, (b) For the proof of Lemma 1 151 a fiber cannot 
have more tlian one intersection point witli purple arcs. 

are pure and a vertex incident to at least two pure arcs pointing toward it is of type (e) or (f), which is 
mixed (see Fig. [3]). Hence, w is a mixed vertex of S{<S), which concludes the proof. □ 

Lemma 14. Given the farthest-polygon Voronoi diagrams of two families Si and S2 of pairwise disjoint 
polygonal sites, the mixed vertices of ^{Si U S2) can be computed in time 0(rilog^ n). 

Computing the mixed vertices in time 0{nlog'^ n) is the most subtle part of the algorithm and we 
postpone the proof of this lemma to Section 14.11 after showing how to compute the purple chains from 
some starting point fLemma llOp . We start with a preliminary lemma which is an important consequence 
of Lemma [S] 

Lemma 15. Let fx be the fiber of point x in a cell C of^{Si) or 5(52). Then, the relative interior of 
fx intersects the purple arcs of ^{S) in at most one point. 

Proof. Assume, for a contradiction, that the relative interior of intersects two purple arcs in two 
distinct points p and q, where q lies on fp (see Fig. |6(b)[ ). We can assume, without loss of generality 
that the purple chain and fx cross at p, because, otherwise, there exists another fiber close to fx (for 
instance, fy for some y close to x) that intersects transversally the purple chains in two points. Now, let 
P be the site containing the feature w associated with C. In C, the purple arcs bound the cell of ^{S) 
belonging to feature w. Hence, there is a point p' on fx sufficiently close to p such that p' G Tl{P) in 
5(5). Moreover, the fiber fpi is, by definition, a subset of the fiber fx. Thus, the fiber fpi contains q 
and thus intersects a purple arc, contradicting the fact that fp' lies in TZ{P) in ^{S) (by Lemma □ 

Lemma 16. Given the farthest-polygon Voronoi diagrams S{Si) and ^{S2) of two families Si and S2 
of pairwise disjoint polygonal sites, the purple chains of^{Si LIS2) can be computed in time 0(rilog^ n). 

Proof. By Lemmas 1121 1131 and 1141 we can compute the vertices at infinity of the purple chains, and a 
superset of size 0{n) of at least one mixed vertex per bounded component. As we have seen in the proof 
of Lemma [T^ these latter vertices are of type (e) or (f) (see Fig. [21); they thus involve only two sites and 
we discard all those that involve two sites of Si or two sites of 52- We thus obtain a set of 0{n) purple 
vertices with at least one such vertex on each purple chain. We consider each such vertex, denoted v, in 
turn, and "trace" (construct) the purple chain that v lies on. 

If V is at infinity, there is only one purple arc incident to it; otherwise, there are two and we consider 
one of them. The bisector supporting the incident purple arc is that of one red and one blue feature 
among the features defining v. Splitting the bisector at v defines two semi-infinite curves incident to v 
and we can determine, using the respective position of the three features defining v, which of these two 
curves supports the considered purple arc; we call this semi-infinite curve a purple half-bisector. 

We then trace the purple chain by following its purple arcs from cell to cell: observe that the other 
endpoint of the considered purple arc incident to v is either at infinity or is the first intersection point 
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(starting from v) between the purple half-bisector and the cell boundaries of 5'(5i) and ^{S2)- We 
compute this point as follows. 

We first locate v in S{Si) and ^{S2)- This can be done in O(logn) time, assuming that we have pre- 
computed a point-location data structure for ^(Si) and 5'(52) in 0(n log n) time (see, for instance, |3).EI 
Note that w is a vertex of S{S) and that it involves two features of one site, say in Si, thus v lies on a 
medial axis arc of induced by these two features. We then also determine the cell of ^{Si), denoted 

Ci, that contains the purple half-bisector in a neighborhood of v (this is a constant size problem). Since 
t; is a purple vertex, its third feature belongs to a site of 52, and v lies in the cell of 5'(52), denoted C2, 
belonging to that feature. Let Wi denote the feature associated with Ci, i = 1,2. 

Then, for i = 1, 2, we sweep the cell Ci with a line orthogonal to the feature Wi if Wi is an edge and 
with a line through the feature Wi if Wi is a corner. If Wi is a corner, the sweep is done clockwise or 
counterclockwise so that the sweep line intersects the half-bisector inside Ci (deciding between clockwise 
or counterclockwise is a constant size problem); the situation is similar when Wi is an edge. 

The two cells Ci and C'2 are swept simultaneously. However, since the two sweeps are not a priori 
performed using the same sweep line, this requires some care. For clarity, we first present each sweep 
independently. 

By Lemma [51 the sweep line always intersects one arc of the upper chain of Ci (or two arcs at their 
common endpoints) and similarly for the lower chain. We first determine the arcs of the upper and lower 
chains that are intersected by the sweep line through v (or about to be intersected if the sweep line goes 
through a vertex of the chain). We also determine the intersection, if any, of these two arcs with the 
purple half-biscctor. When the sweep line reaches an endpoint of one of the two arcs that are being 
swept, we determine the intersection (if any) between the new arc and the purple half-bisector. When 
the sweep line reaches the first of the computed intersection points between the purple half-bisector and 
the boundary of the cell, we report this intersection point and terminate the sweep of Ci . 

If the two sweeps were performed independently, we could report the first point where the purple 
half-bisector exits one of the cells Ci or C2, and continues the tracing in a neighboring cell of either 
^(iSi) or 5^(52), along a new purple half-bisector. This however would not yield the claimed complexity 
because, roughly speaking, if a sequence of purple arcs enter and exit 9(n) cells Ci, C(, C", . • • while 
remaining inside a cell C2 of complexity &{n), the cell C2 would be swept many times, possibly leading 
to a complexity of Q{n\ogn) per arc, and a total of <d{n^ \ogn). We thus perform the two sweeps 
simultaneously, as follows. 

Note first that, by Lemma |151 during the sweep of C'i, the point of intersection, in Ci, between the 
sweep line and the purple half-bisector moves monotonically along the purple half-bisector. We can thus 
parameterize the sweeps of Ci and C2 by a point x moving monotonically on the purple half-bisector 
away from v. The point x define two sweep lines (the lines through x and through Wi or orthogonal to 
Wi depending on the nature of Wi), and the events are those of the sweeps of Ci and C2. This sweep 
ends when the purple half-bisector leaves Ci or C2 . Then, the tracing of the purple chain continues in a 
neighboring cell of either ^{Si) or ^{S2), along a new purple half-bisector. We stop tracing the purple 
chain when we reach a vertex at infinity or the vertex v we started from. 

We then consider a new starting point v; note that we can easily check in 0(log n) time whether it 
has already been computed (while tracing some purple chains) by maintaining the list of the already 
computed vertices on the purple chains, ordered lexicographically by their features. 

We now analyze the complexity of the algorithm. Consider first the initialization of every sweep, 
that is the determination of the arcs of Ci and C2 that are intersected by the two sweep lines through 
V. There are 0{n) such initialization steps to perform, since ^(S) have size 0{n) by Theorem [51 and 
each step can be done using binary search in O(logn) time, after preprocessing all the cells of ^{Si) and 
5^(52); the preprocessing (which simply is storing the ordered vertices of the upper and lower chains of 
each cell in arrays) can be done in time linear in the total size of the cells, which is 0{n), by Theorem [3 

We finally analyze the complexity of the rest of the algorithm by applying a simple charging scheme. 
Note first that intersecting the purple half-bisector with an arc of a cell takes constant time. We charge 
the cost of computing these intersections to either the purple vertices or to the arcs of Ci, as follows. The 
intersections with each of the arcs that are swept at the beginning and the end of the sweep are charged 
to the corresponding endpoint of the purple arc. Every purple vertex is thus charged at most eight times 

^In fact, we will see in Section |4. II that the combinatorial description of v and all the information regarding its location 
in 5{<Si) and S{S2) are already available as a by-product of the computation of the mixed vertices. This makes the location 
procedure described here not strictly necessary. 
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because each of the two incident purple edges is intersected with one arc of each of the lower and upper 
chains of each of the two cells Ci and C2 . The intersections with each of the other arcs of Ci are charged 
to the arc in question. Every such arc is charged at most once per cell, and thus at most twice in total; 
indeed, such an arc of Ci is swept entirely during the sweep of Ci and thus, by Lemma [T51 it will not be 
swept again during another sweep of cell Ci (when treating another purple half-bisector). Since ^{Si), 
5'(iS2)j and 3^(5) have size 0{n) (Thcorcm|n|), all the purple arcs can be computed in 0{nlogn) time, in 
total, once the set of starting points are known, and thus in 0(nlog^ n) time by Lemmas [T^ and [HI □ 

Lemma 17. Given the farthest-polygon Voronoi diagrams of two families Si and S2 of pairwise disjoint 
polygonal sites, the farthest-polygon Voronoi diagrams of Si U S2 can he computed in time 0{n\o^ n). 

Proof. During the computation of the purple chains (see the proof of Lemma I16p , we can split the arcs 
of 5^(5i) and 5'(iS2) at every new purple vertex that is computed. Then, merging ^{Si) and ^(52) can 
trivially be done in linear time since, as we mentioned before, the diagram iJ(iS) consists of those portions 
of S'(5i) lying in the red regions of 5^(5), and those portions of ^{82) lying in the blue regions of S'(iS), 
see Fig. [S] □ 

4.1 Computing the mixed vertices 

We prove here Lemma |14I stating that, given the farthest-polygon Voronoi diagrams of two families iSi 
and S2 of pairwise disjoint polygonal sites, the mixed vertices of ^{Si U1S2) can be computed in time 
0(nlog^ n). 

We start by computing the randomized point-location data structure of Mulmuley [T7] (see also [3 
Chapter 6]) for the two given Voronoi diagrams t?(5i) and ^{82) ■ (We chose this randomized algorithm 
for its simplicity; randomization can however be avoided as we explain at the end of this section.) This 
data structure only needs two primitive operations: (i) for a given point p in the plane, determine 
whether the query point lies left or right of p, and (ii) for an x- monotone line segment or parabolic arc 7, 
determine whether the query point lies above or below 7. Both cases can be summarized as follows: 
Given a comparator 7, determine on which side of 7 the query point lies. The comparator can be cither 
a line or a parabolic arc. 

We compute the mixed vertices lying on each tree T of each medial axis 9Jl(P) separately. The 
intersection T H 7^(-P) is a connected subtree by Lemma We can locate the internal vertices of 

this subtree easily, by performing a point location operation for each vertex u of T in 5^(<Si) and 5'(iS2), 
deducing which site is farthest from v and checking if the farthest site from v is P. Let / be the set of 
vertices of T that lie in TZ{P). We now need to consider two cases. 

4.1.1 When / is not empty 

If / is non-empty, then every arc a of T incident to one vertex in / and one vertex not in / must contain 
exactly one mixed vertex s* by Lemma [2 To locate this vertex s*, we use parametric search along the 
arc a [14] , albeit in a restricted way, as we do not need to use parallel computation: The idea is to 
execute two point location queries in ^{Si) and ^{82) using s* as the query point. Each query executes 
a sequence of primitive operations, where we compare the (unknown) location of s* with a comparator 
7 (a line or a parabolic arc). This primitive operation can be implemented by intersecting a with 7, 
resulting in a set of at most four points. In O(logn) time, we can test for each of these points whether 
it lies in TZ{P). This tells us between which of these points the unknown mixed vertex s* lies, and we 
can answer the primitive operation. 

It follows that wc can execute the two point location queries on s* in time 0(log^ n), and we obtain 
one cell and one arc of ^{Si) and ^{S2) containing s* (for instance if P G iSi then s* lies on an arc of 
d{Si) and in a cell of ^{S2))- The mixed vertex s* lies at equal distance of the three features to which 
the arc and cell belong (two of which arc features of P) . 

4.1.2 When / is empty 

It remains to consider the case where / is empty, that is, no vertex of T lies in TZ{P). Nevertheless, the 
region TZ{P) may intersect a single arc a of T. In this case. Lemma [3] tells us that there are two mixed 
vertices on a that we need to find. We first need to identify the arcs of T where this could happen. 



13 



Figure 7: (a) pq is a sub-arc of 9Jl(P). p'q' is the intersection of pq with TZ{P). We also have p'q' ~ 
97l(P) n7^(P). (b) The cyhndcr €{p, q) of the pair (p, g) in (a), (c) Cyhndcrs and £((7, r) illustrate 

the two other possible geometric shapes of a cylinder, whose boundary alternate between a part of a 
feature of P and a circular arc. 

Let p and g be two points on the same arc a of *It(P). We define the cylinder €{p,q) of the pair 
{p, q) as 

where the union is taken over all points x on the arc a between p and q. Since p and q belong to the 
same arc a of 9JT(P), which is induced by either two edges, two vertices, or one edge and one vertex of P, 
cylinders may have only three possible shapes, illustrated in Fig. [Tl^b, c). We define a condition G{p,q) 
as follows: Let Q be a site farthest from p, and let w be a feature of Q closest to p. Then G{p, q) is 
true if u; C <^{p,q) or if Q intersects Dp{q). Note that G{p,q) could possibly depend on the choice of Q 
and w when the farthest site or the closest feature is not unique. We will show below that this choice 
does not matter for the correctness of our algorithm. 
We can now prove the following two lemmas: 

Lemma 18. Let p,q be points on the same arc a of^Xll{P), such that neither p nor q lie in TZ{P). If a 
intersects 7?.(P) between p and q, then G{p, q) and G{q,p) both hold. 

Proof. Since the statement is symmetric, we only need to show G(p, q). Let Q be a site farthest from p, 
and let x be a point between p and g on a that lies in TZ(P). This implies that Q intersects Ds{x) = 
Dp{x). Since Q does not intersect Dp{p), we deduce that Q intersects €{p,q). If Q does not lie entirely 
in q) then, since the sites P and Q are disjoint, Q must cross the common boundary of €{p, q) and 
Dp{q) (see Fig. [Tl^b, c)). Thus Q intersects Dp{q) and indeed G{p,q) holds. □ 

Lemma 19. Let p, q be points on the same arc a of a tree T of 9Jl(P) that admits no vertex in TZ{P). 
If neither p nor q lie in TZ{P) and both G{p, q) and G{q, p) hold, then all points in T D 7?.(P) lie on a. 

Proof. Assume, for a contradiction, that there exists a point x on T in TZ{P) not between p and q on a; 
assume, without loss of generality, that p lies on the path joining x and q in T. Since G(p, q) holds, there 
is a farthest site Q ^ P from p such that Q intersects q). Therefore there exists a point y between p 
and g on a such that y is closer to Q than to P. Since x S 7?,(P) by assumption, x is closer to Q than 
to P. We have shown that x and y are closer to Q than to P but the point p, which is between x and y 
on T, is closer to P than to Q. This contradicts Lemma [1] and concludes the proof. □ 

Let us call an arc a connecting vertices p and q of T a candidate arc if G{p, q) and G{q,p) both hold. 
Lemma [T51 implies immediately that if there are two candidate arcs, then T Ci 'R-iP) is empty, and there 
are no mixed vertices on T. 
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Since we have point-location data structures for 5'(<Si) and S'(iS2), we can test the condition G{p,q) 
in time O(logn) for a given arc a in 9Jl(P) and two points p,q G a. This aUows to identify aU candidate 
arcs in 0{mlogn) time, where m is the complexity of T. If there are zero or more than one candidate 
arcs, we can stop immediately, as there are no mixed vertices on T. 

It remains to consider the case where there is a single candidate arc a in T. We again apply parametric 
search, using an unknown point s* in a H 'R-{P) as the query point. During the point-location query, we 
maintain an interval on a that must contain s* (if s* exists at all). To implement a primitive query, we 
must determine the location of s* with respect to a comparator 7. If the current interval pq on a does not 
intersect 7, we can proceed immediately, otherwise we get a sequence of points p = xq, xi, X2, . . . , ^ q 
on a, where 2 ^ s ^ 5 (since there are at most four intersection points xi, X2, ■ ■ ■ , Xg^i between two 
curves of degree at most two). We first test if some Xi £ 7?.(P) in C>(logn) time. If so, we abort the 
process, and use the method discussed above to find the mixed vertices on the arc between p and Xi 
and between Xi and q. If no Xi lies in TZ{P), we test the conditions G{xi,Xi+i) and G{xi+i,Xi) for each 
consecutive pair. If the conditions hold for no pair or for more than one pair, we can stop immediately, 
as there cannot be a point of 7?.(P) on a. If there is exactly one pair, we have found the location of s* 
with respect to the comparator 7, and we continue the point location query. 

If both point location queries on s* in 5^(5i) and '^(82) terminate without encountering a point 
in Tl{P), the current arc pq of a lies entirely on an arc of i?(5i) and in a cell of 5'(iS2) (assuming that 
P G Si). Moreover, we get this arc and cell from the two point location queries. The two mixed vertices 
s* on the arc pq are then both defined by the same three features that define this arc and ccll|f| and they 
can be computed in constant time. 

4.1.3 Avoiding randomization 

Finally, we argue that, instead of using a randomized point-location data structure in the above algorithrn, 
we can use any other data structure, such as the one of Edclsbrunncr et al. 0, as long as all predicatct|j 
used in the associated point-location algorithm arc answered by evaluating the signs of polynomial 
expressions of bounded degree in the input data. If one predicate is answered by evaluating the sign 
of several polynomial expressions, it can be split into several predicates, each of which corresponds to 
exactly one polynomial expression. Then, a predicate corresponds to a polynomial expression of bounded 
degree that depends on the x and y coordinates of the query point and a set of other parameters (for 
instance, the coordinates of a point, or the coefhcients of an implicit equation of a curve, against which 
the query point is tested); this expression, seen as a polynomial in x and y, defines a curve, 7, of bounded 
degree. Recall now that we perform a point location query using an unknown query point s* that lies 
on a straight or parabolic arc. The curve 7 associated with a polynomial predicate splits this arc into a 
bounded number of pieces along which the sign of the polynomial is constant. To answer the query, it 
suffices to use the method described above on the boundary points of these pieces. 

5 Concluding remarks 

We have considered, in this paper, farthest-site Voronoi diagrams of k disjoint connected polygonal 
sites in general position and of total complexity n. In particular, we proved that such diagrams have 
complexity 0{n) and that they can be computed in 0(ri, log^ n) time. 

We have seen that Voronoi regions can consist of several unbounded components. However, since the 
pattern PQPQ cannot appear at infinity, it is always possible to connect the components of one Voronoi 
region by drawing non-crossing connections "at infinity," and so we can think about Voronoi regions as 
being connected at infinity. This curious property can perhaps be better understood by studying the 
same problem on the sphere. Here the resulting structure is simpler: The bisector of two polygons is 
a single closed curve, and the family of bisectors of a fixed polygon P with the other polygons forms 
a collection of pseudo-circles. (The closest-site Voronoi diagram of three disjoint polygons cannot have 
three vertices.) The farthest-site Voronoi region of P is the intersection of the pseudo-disks that do not 

''Three features may define two mixed vertices, for instance in tlie simple special case where the sites consist of a 
V-shaped polygonal site and a point. 

^For instance, in the case of Mulmuley's point- location data structure 1171 . the predicates (i) and (ii) mentioned at the 
beginning of Section |4. II 
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contain it, and is thus either empty or simply connected [16| . The 0{k) bound on the number of pure 
vertices is then a simple consequence of the planarity of this diagram. (With some care, an alternate 
proof of Lemma [5] based on this pseudo-disk property could be given.) 

Farthest-site Voronoi diagrams are related to the function 
fp : X ^ argmax5g5 (minj,g5 (i(a;, y)) which returns the farthest site to a query point. The standard 
closest-site Voronoi diagram corresponds to the function a; ^ argminsg^ (min^gs d{x, y)) which returns 
the closest site to a query point. Voronoi diagrams induced by other similar functions have also been 
considered in the literature. In particular, the diagram obtained by considering the "dual" function of 
fp, that is a; I— >■ argminggg (maxj,g5 (i(a;, y)) is the so-called Hausdorff Voronoi diagram; see [18j and 
references therein. 

A Hausdorff diagram is typically defined for a collection of sets of points in convex position because 
the maximum distance from a point to a polygonal site is realized at a vertex of the polygon's convex 
hull. Papadopoulou showed that the size of the Hausdorff Voronoi diagram is 8(71-1- M), where n is 
the number of points in the collection and M is the number of so-called "crucial supporting segments" 
between pairs of "crossing sets" (a pair of sets is crossing if the convex hull boundary of their union 
admits more than two "supporting segments", that is, segments joining the convex hulls of each set; such 
a segment is said crucial if it is enclosed in the minimum enclosing circle of each set.) 

A fast parallel algorithm was obtained by Dehne et al. for the special case when Af = ^7^. They 
gave an 0{{n\o^n)/p) time parallel algorithm for the diagram construction on p processors and a 
0(n log'* n) time sequential algorithm. Their algorithm is similar to ours and differs mainly is the way 
the purple chains are constructed. Indeed, the arcs of a Hausdorff Voronoi diagram are (straight) line 
segments; this permits the use of an ad hoc data structure for finding the "mixed" vertices. In contrast, 
arcs in a farthest-polygon Voronoi diagram can be curved and we offer a technique for computing the 
purple chains using parametric search, which is more efficient and more general. This generality has 
already found another application in the construction of farthest-site Voronoi diagrams for the geodesic 
distance [S] . We also believe that our technique could be applied in the sequential algorithm of Dehne et 
al. to improve its time complexity to O(nlog^n). 

On the other hand, we do not provide a parallel algorithm for computing farthest-polygon Voronoi 
diagrams. It would be of interest to study the feasibility of applying the parallel techniques of Dehne et 
al. to the farthest-site Voronoi diagram computation. 
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